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Abstract: Linear rate equations are used to describe the cascading decay of an initial 
heavy cluster into fragments. Using a procedure inspired by the similar, but continuous 
case of jet fragmentation in QCD, this discretized process may be analyzed into eigen- 
modes, corresponding to moments of the distribution of multiplicities. The orders of these 
moments are usually noninteger numbers. The resulting analysis can be made time inde- 
pendent and is applicable to various phenomenological multifragmentation processes, in 
which case it leads to new approximate finite-size scaling relations for the spectrum of 
fragments. 
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In this work we consider binary fragmentation processes where any fragment with 
mass number k breaks into fragments with mass numbers j and k — j, j = 1, 2. ..A; — 1 

m 

with a probability Wjk per unit of time. It is assumed that Wjk is time independent. 
By definition, Wjk is symmetric if j is replaced by k — j, naturally. Let Nj(t) be the 
multiplicity of fragment j at time i in a process initiated from the decay of a cluster A, 

r- 1 | 

namely iVj(O) = 5ja- The model under study is described by the following set of linear, 



first order differential equations, 



—L = - CjNj + J2 w jk N k , j = 1, ...A, (1) 
k=j+l 

with 

With components Nj, j = 1, ...A, for a column vector |jV >, the system, Eqs.(l), obviously 
boils down to d\J\f > /dt = W|jV > with a triangular matrix W. For the sake of clarity we 
show here W when A = 4, 



/0 W 12 Wis t"l4 \ 

—W 12 /2 W 2 3 = U>13 W 2 A 

-(wi3+W 2 3)/2 ^34 = ^14 

\0 -(^i4 + W 2 4 + M'34)/2/ 



(3) 
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The general solution of Eqs.(l) is obviously a sum of exponentials whose rates of decay 
in time are the diagonal matrix elements — Cj , the trivial eigenvalues of the triangular W. 
(We notice that any increase of the dimension A leaves intact the preexisting eigenvalues 
and only adds new ones.) 

This matrix W has a remarquable property, namely a fixed left (row-like) eigenstate 
All, whose components Alij = j, j = 1, ...A, do not depend on the wg^s. This comes from 
the symmetries of W demanded by the conservation of the total mass Mi = ^2f = ijNj, 
with dMi/dt =< Ali|W|A > . The corresponding eigenvalue is, naturally, —c\ = 0. 
Moreover, it is clear that the other eigenvalues, — ca, A = 2, ...A, induce a triangular 
matrix of left" 1 " eigenstates M\, namely the components M\j vanish when j < A. 

The higher A, the more the corresponding (bra-like) eigenstate is only probing heavy 
fragments in the multiplicity distribution A/". This is reminiscent of a hierarchy of mo- 
ments M q = ^f = ij q Nj, when the exponent q increases. Even though moments do not 
make a strictly triangular rearrangement of the information contained in A/", they rep- 
resent a natural continuation of the first eigenvector All, and this letter will show that 
there is a practical connection between moments and eigenvectors. Indeed, besides trian- 
gularity, there is another argument indicating the interest of moments in the solution of 
rate equations like Eqs.(l). In field theory models of jet fragmentation at high energy, 
similar evolution equations appear as solutions of perturbative QCD; they are continu- 
ous, since the quantity which fragments is the energy-momentum of quarks and gluons, 
before their transformation into the observed hadrons. These equations are well known 
as Gribov-Lipatov-Altarelli-Parisi (GLAP) equations^. As a function of the fraction x 
of energy-momentum, the equations for the fragmentation function N(x) can be exactly 
diagonalized by moments M q = J dxx q N{x). In this case, all values of q > 1 are admitted, 
as a consequence of the continuous character of the equations. Our task, in this letter, 
is to examine the effect of discretization on this field-theoretical result. Note that such a 
discretization have been introduced long ago^, but only for a numerical approximation of 
the continuous equations. 
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In the following, we will consider "triangularly redefined" moments 



A 



(4) 



where the set of exponents q(X), A = 2, ...A, may contain nonintegers. The time derivative 
of such a redefined moment is 



where we have used an interchange of indices j and k in the double summation, and also 
the convention that diagonal rates Wkk are identically vanishing. The double summation 
is sketched on Fig.l, displaying the weight matrix Wjk- 

It turns out that, in practice, there exists special values q(X) for which the coefficients 
d(X, k) happen to be, at least approximately, independent of k when k is large enough. 
Then, they can be factored out in formula (5) and they appear as eigenvalues for the corre- 
sponding eigenvector M\. Moreover, simultaneously, since we know the exact eigenvalues 
(see Eq.(2)), one gets d(X, k) ~ — c\. This will be shown both analytically in a "continu- 
ous" limit (i.e. large matrices) and numerically for various models with finite sizes A. In 
fact one is led to consider first the discrete models which admit the field-theoretical type 
of equations in the continuous limit - let us call them the scale-invariant case, since no 
dependence on the matrix size appears explicitely, and then the more general situation. 

1. The scale-invariant case 

When A is large, A finite, and k large but smaller than A, then the ratio x = j/k 
can be considered as a continuous label, < x < 1, in Eq.(6). Moreover, with positive 
values of q, an extension of the first summation in Eq.(6), $^j=A' into a summation Y^=i 
brings a weak contribution from < x < X/k. From QCD, where x corresponds to the 
fragmentation of momentum and a scale-invariant property of the transition weights is 





with 
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valid, one may consider a large class of models setting 



Wjk = <fi(j/k)/k = <fi(x)/k, 



(7) 



where (p is any suitable function of the scaling variable x. More precisely, because of 
the symmetry necessary for Wjk, a large class of legitimate models correspond to Wjk = 



[f(j/k) + f(l - j/k)]/(2k) = [f{x) + f(l - x)]/(2k). Hence, for large values of k, both 



where dx replaces 1/k. It will be noticed that d(X) does not depend any more on k. It still 
depends on A via the exponent <?(A), naturally. Note also that the "splitting" function / 
can be general, even with some singularity at both ends of the integration domain, provided 
the integral itself converges. 

There may also be a continuous limit for c\ if A becomes large. Indeed, according to 
Eq.(2), c\ — > Jq dxip(x)/2, if 1/A amounts to dx and if this integral converges. In such a 
case, the spectrum of W accumulates into a quasi degeneracy. However, it is important to 
realize that the convergence of the c^'s is not required. On the contrary, it is quite possible 
that the limiting continuous model does not exist, leading to an infinite hierarchy of q(X). 
We shall meet such cases later on. Moreover, for low values of the label A, this continuous 
limit is not in order. 

We now notice that (f(x) is a semi-positive definite function since Wjk, a transition 
rate, cannot become a negative number. Hence d(X) is a monotonically decreasing function 
of q. It vanishes for q = 1, as expected from the conservation of M x . According to Eq.(5), 
the time evolution of a triangular moment M\ becomes very simple if d(X) can be identified 
with the eigenvalue, — c\. Hence, for each A, we consider the exponent q(X) which is the 
unique solution of the consistency equation 



The discrete set of solutions q\ of this equation, when the integer label A runs from 
1 to A, define the "eigenmoments" of the theory, namely those moments whose time 
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summations in Eq.(6) amount to the discretization of an integral 




(8) 





evolution is (almost) proportional to just one exponential exp(— c\t), rather than a mixture 
of such exponentials. It will be noticed that, since the sequence of coefficients c\ increases 
monotonically, the sequence of solutions q(X) is also a monotonically increasing sequence, 
starting from q(l) = 1 with c\ = 0. 

Set temporarily f(x) = 1/x 13 ', with (3=1. This case is reminiscent of the QCD evolu- 
tion equations for gluons whose kernel contains the same singularity at small x W. Then 
Wjk = [1/j + l/(k — j)]/2, and one finds easily that C2 = 1/2, and that the corresponding 
solution of Eq.(9) is q(2) = 2. One also finds that C3 = 3/4, and that the corresponding 
solution of Eq.(9) is g(3) = 3. More generally, one finds for Eq.(9) the solution q(X) = X. 
This definitely suggests that integer moments form an infinite sequence of eigenmoments 
for that choice of /, f(x) = 1/x. In agreement with this hint, ones finds easily from Eq.(6) 
that the sequence of coefficients d(2, 20) = -0.475, d(2, 21) = -0.476,...d(2, 49) = -0.4898, 
d(2, 50) = —0.4900,... converges towards — ci. Just to give another example, the sequence 
d(5,20) = -0.980, d(5,21) = -0.983,. ..d(5, 49) = -1.016, d(5,50) = -1.017,... converges 
towards — C5 = —1.042. And so on for all the moments M\, which thus generate excellent 
approximations to eigenvectors when the exponents q(X) are just integers. 

This argument is independent from the normalization of /, since Eq.(9) is homo- 
geneous with respect to trivial multiplications of / by an overall constant. It may be 
interesting to note that the diverging sequence of values for q(X) might be related to the 
divergence of the eigenvalue sequence at infinite A. 

A similar result can be observed if (3 7^ 1, but now the solutions of Eq.(9) do not 
correspond to integer exponents. For instance, with (3 = —0.5, one finds C2 = 0.177, 
q{2) ~ 3.13, c 3 = 0.232, q(3) ~ 5.26, c 4 = 0.259, q(A) ~ 7.4, c 5 = 0.275, q(5) ~ 9.6,... 
The convergence of the coefficients d towards the eigenvalues is still surprisingly good. For 
instance d(2, 20) = -0.176, and rf(5, 20) = -0.275. 

This discussion (for the scale invariant case) is illustrated on Figs. 2. Here, rather 
than asking whether a moment may behave like an eigenvector, we consider the reverse 
question: given an eigenvector M.\, does it happen that the components A4\j induce 
effective moments M q ^l Namely, is there an exponent q(X) compatible with M.\j oc 
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j q ^ X \ j > ^? I n Figs. 2, for instance, the components of the second (A = 2) and fourth 



(A = 4) eigenvectors are displayed as functions of the fragment size. It will be stressed 
that they are almost linear in a Log-Log plot in a large interval starting from the maximal 
choosen value A = 30. These figures, Figs. 2, show the structure of moments depending on 
the parameter (3. 

In table I we show the comparison between the actual eigenvalues c\, (for A = 2 and 4) 
with those obtained through Eq.(9) after the determination of the effective values q(X) from 
Figs. (2). The agreement is pretty good, except perhaps for (3 ~ 2, which is at the borderline 
of convergence of the integral in Eq.(9). It can be thus claimed that various choices of (p 
in Eq.(7) make it possible to find eigenmoments via the continuous limit and solutions of 
Eq.(9). It will be noticed that this continuous limit is closely linked to the denominator 
k in Eq.(7), since this denominator induces the needed measure dx, independently of the 
overall scale given by A. 

2. Scale- dependent cases 

In a more general situation, Wjk is not compatible with Eq.(7). This will in general 
lead to an introduction of the overall scale A in the problem, and corresponds to cases 
where the scale invariance of the weights is not preserved. Let us illustrate this by the 
following instance: 



where a may be different from 1. 

To be specific, but as an example of more general value, we display on Fig. 3 the 
Log-Log plot M. \j versus j, for A = 2, A = 50, f3 = and various values of a. From this 
figure, one realizes on this simple example that the diagonalization by eigenmoments is 
obtained for a < 1, while for a > 1 there is a clear distorsion of eigenvectors with respect 
to moments. It is interesting to interpret this phenomenon analytically by inspecting the 
modification which occurs with the choice of Eq.(10) for the weights. Instead of Eqs.(5-6), 
one finds the following, 



[V/k)-e + (i-j/k)-i>]/(2kr), 



(10) 



dM x 



dt 




(11) 
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where w corresponds to rescaled weights with a = 1. Following the previous discussion, 
based on the existence of a continuous limit <i(A, k) ~ <i(A), one is led to the approximate 
consistency equation: 

-d{\) ~ A a_1 c(A), (12) 

where the renormalisation factor takes care of the initial values of the moments. 

The occurrence of the A— dependent factor is the signal of the lack of scale invariance 
of the fragmentation dynamics when a ^ 1. Note that Eq.(12) uses the assumption that 
the eigenvalues do not change substantially between q and q — 1 + a. It can only be an 
approximation. 

In Table II, we display the different values obtained for the q(2) for various values of 
a, obtained by the fitted slopes at the origin for the curves obtained in Fig. 3. We compare 
them to those obtained from Eq.(12) when the input are the actual exact eigenvalues c^- 
The agreement is here also quite satisfactory, except in the region when a > 1. One might 
associate this phenomenon to the well-known! 3 ! fact that a shattering transition takes place 
at finite time, the conservation of mass being broken in the continuum limit. A special 
study of this case is in order for the future. 

The class of models which can be analyzed by eigenmoments is thus larger than the 
class described by Eq.(7). It must be noted, however, that the "eigenorders" q(X) are not 
universal, but clearly model-dependent. 

As an application of the properties of eigenmoments, let us consider the problem of 
3-dimensional bond percolation in a finite-size square lattice. This model seems to give 
a successful description of nuclear multifragmentation, when a heavy ion receives enough 
excitation energy to form a highly unstable state and decays into several fragments^ . The 
statistics of fragment numbers and sizes seem to follow predictions of a percolation model 
in which each lattice site is populated by a nucleon, and the percolation parameter p, 
namely the survival probability for bonds, varies between and 1. There is no obvious 
time scale in this model, hence a comparison of its predictions with those of linear rate 
equations models requires the use of eigenmoments in order to obtain an intrinsic time 
scale from the evolution of such eigenmoments. 
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For this purpose, we remark that, if they are identified as eigenmoments, the M 9 's are 
linked by linear relations in Log-Log plots, and their explicit time dependence disappears. 
We are thus led to display in the same way the moments obtained from the percolation 
model, choosing for instance M2 for reference, see Fig. 4. Different moments are displayed 
(with q = l, 1.5, 2, 3, 4, 5) and show the interesting feature of a quasi-linear dependence for 
the values q = 3, 4, 5, given the fact that for q = 1 (mass conservation) and q = 2 (reference 
scale) the linear dependence is fixed. It is clear from this figure that the quasi-linear form 
is obtained between p = 1 and p = p c , where p c is the critical value above which, in the 
continuous limit, an infinite percolation cluster is formed. Indeed, the figure shows the 
dominant contribution of the cluster of largest mass to the averaged moments. This largest 
cluster is, for finite size problems, the representative of the infinite cluster when p > p c . 

Notice that the moments implied by the rate equations are the full moments, including 
the largest fragment, while in usual analyses of percolation models^, scaling properties are 
investigated with moments modified by the subtraction of the largest cluster. Moreover, in 
such traditional analyses of percolation, the reference time scale is generally given by the 
moment Mq or a similar variable related to the multiplicity of fragments. The comparison 
and compatibility of our approach with such analyses is an open problem of some interest ^ . 

In conclusion, from this first study on the percolation model, we obtain a hint that 
linear rate equations could provide a time dependent description of multifragmentation. 
But it is difficult at this stage to obtain informations on the set of eigenorders q(X) which 
could be associated with percolation. The existence of scaling relations between moments 
can be proven in the vicinity of p c for percolation through finite-size scaling^ 4 !. Our result 
is compatible with this and, furthermore, involves the whole region 1 > p > p c . In the 
representation provided by rate equations, we have obtained scaling relations valid for 
short time scales, while the previous general results, see for instance RefJ 7 ', involve long 
time scales only. An open problem is to connect both analyses for a general system of 
equations. 

Acknowledgments R.P. thanks Xavier Campi for his patient explanations on nuclear 
multifragmentation and for providing the authors with the suitable percolation program. 
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^Footnote The right-hand-side, ket-like eigenstates define also a triangular matrix, 
apparently unrelated to the matrix of bra eigenstates, except for trivial biorthogonality 
relations. Up to now we have been unable to find a practical use of these ket-like eigen- 
vectors. 
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Figure captions 
Fig.l: Discrete rate equations: Double Summation 

The double summation range for truncated moments M\ is represented by the hatched 
triangle. White dots: diagonal weigths Wjj are zero. Black dots: non-diagonal weigths 
Wjk- The interchange of indices j and k in the description of the hatched triangle leads to 
Eqs.(5-6). 

Fig. 2: Eigenvectors: scale- invariant case (a = 1) 

The components of the rate equation eigenvectors M\j are displayed as functions of j in a 
Log-Log plot. The chosen weights correspond to Eq.(lO) with a = 1, and different values 
of (3. The curves correspond to a smooth interpolation (dark line), resp. extrapolation 
(dashed line), of the exact eigencomponents for a system of size A = 30. Fig 2-a: A = 2 ; 
Fig 2-b: A = 4. The curves are used for the determination of the effective values q(2), (/(4), 
see Table I. 

Fig. 3: Eigenvectors: scale-dependent cases {(3 = 0) 

Same as Figure 2-a but for weights following Eq.(10) with (3 — and different values of a. 

The curves give the eigenorders listed in Table II. 

Fig. 4 Percolation analyzed with the "M 2 time scale" 

Relative strengths of moments M q , q = 0, 1, 1.5, 2, 3, 4, 5, as functions of M 2 in a Log-Log 
plot. Data taken from 3-d bond percolation on a 6 * 6 * 6 lattice. The corresponding values 
of the bond survival probability p are shown on the horizontal axis. Its critical value is 
p c = .25. Full lines: moments. Dashed lines: contributions of the largest cluster. Dashed- 
dotted line: the reference moment M 2 . Notice that a linear behaviour is approximately 
obtained for < p < p c and q = 3,4, 5. 

Table Captions 

Table I: For a = 1, A = 2, 4 and different values of f3, the effective exponents q(X) from 
Figs.2-a,b and the corresponding values for — d(X), see Eq.(9). The latter are compared 
with the exact eigenvalues c\. 

Table II: For (3 = 0, A = 2 and different values of a, comparison of the measured 
effective eigenorders (obtained from Fig. 3) with those predicted from Eq.(9). 
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p 


</(2) 


-d(2) 


C2 


g(4) 


-d(4) 


c 4 


-2 


2.4 


.061 


.062 


6.1 


.11 


.11 


-1.5 


2.7 


.09 


.09 


6.5 


.14 


.14 


-.5 


3.1 


.18 


.18 


7.4 


.26 


.26 





3 


.25 


.25 


7.1 


.38 


.38 


.5 


2.6 


.36 


.35 


5.8 


.57 


.57 


1 


2 


.51 


.50 


4.2 


.94 


.92 


1.5 


1.5 


.79 


.71 


2.8 


1.8 


1.6 


2 


1.1 


.85 


1 


1.7 


4.6 


2.7 



Table I 



a 


9(2) 
measured 


predicted 





1.03 


1.04 


.5 


1.17 


1.22 


.6 


1.26 


1.33 


.8 


1.64 


1.7 


1 


3.0 


3 


1.1 


5.7 


5.2 


1.2 


20 


14 



Table II 
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